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ABSTRACT 

The problem of the formation of the Moon is still not explained satisfactorily. While it is a 
generally accepted scenario that the last giant impact on Earth between some 50 to 100 million 
years after the starting of the formation of the terrestrial planets formed our natural satellite, 
there are still many open questions like the isotopic composition which is identical for these 
two bodies. In our investigation we will not deal with these problems of chemical composition 
but rather undertake a purely dynamical study to find out the probability of a Mars-sized body 
to collide with the Earth shortly after the formation of the Earth-like planets. For that we 
assume an additional massive body between Venus and Earth, respectively Earth and Mars 
which formed there at the same time as the other terrestrial planets. We have undertaken 
massive n-body integrations of such a planetary system with 4 inner planets (we excluded 
Mercury but assumed one additional body as mentioned before) for up to tens of millions 
of years. Our results led to a statistical estimation of the collision velocities as well as the 
collision angles which will then serve as the basis of further investigation with detailed SPH 
computations. We find a most probable origin of the Earth impactor at a semi-major axis of 
approx. 1.16AU. 
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1 INTRODUCTION 


An assumed giant impact of an additional Mars-sized object 
(Theia) onto the Earth could have led to the formation of the Moon 
after the planets already had their actual mass and no more gas 
was present in the Solar System. Many recent publicati ons dea l 
with this topic, e.g., 


ar System. Many recent publications deal 
Asphaug 120141). Njdayma^fcStevensonJ 


d20j^)^_[(^uaiTes_^^Liss a uerJ__^ d2015h . |jacobson_eta]_ 
Ijacobson & Morbidelli I d2014h [=JM1, Izidoro_etal 


the first ideas dev eloped by [Hartmann & Davis 



ICameron & Ward! j 1 976l ). and later Canup & A sphau 
Detail e d collis i on sc e naria were s tudied e.g., by ICameron I 
dl997l) . ICanup I d2004l) . ICanup I (2008), and ICanup I d20 1 3h where 
the collision was modelled with the aid of so phisticated SPH 
codes. In a most recent article iKaib & Cowan 1 d2015l) [=KC] 
the authors concentrate on the feeding zone of the planet to 
form with respect to the planet’s volatile inventory and isotopic 
composition. Because of the highly random outcome they ask 
the question of how deterministic the outcome of the planetary 
formation is. In fact the correspondence of the results of the 
different model computations of n-body codes is very small. Most 
of these modelisations have been undertaken to understand the 
architecture of our Solar System, which results only from a subset 
of the chosen initial conditions. In KC the authors estimated the 
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likelihood that the mass of Theia could be equal to the mass of 
the Earth, but the probability is rather low. Their results coincide 
well with JM who estimated the collision probability of bodies 
with comparable masses as being low. According to these results 
we have fixed the mass of the additional planet (the ‘p r °j ect ile 
planet’) to m Mars for our computations. Other investigations aimed 
for high velocity encounters e.g. the one by Cuk et all d2012h 
who assumed high velocity collisions for smaller masses of Theia 
(0.025 m Earth < m Theia < 0.05 /n Eart h), but the results of KC show 
that such an event may not be very probable because a spin rate 
of the Earth of the order of 2 hours can only be achieved by big 
impactors - and these events_are rare. Because of the same reason 
the scenario proposed by [Reufer et al I d2012h where they look for 
a steeper collision angle is not very probable. KC undertake 150 
different simulations with 3 different underlying models: a first 
model with Jupiter and Saturn on circular orbits, a second one with 
initially small eccentriciti es of Jup iter and Saturn, and a third one 
according to the model of lHansenl d2009h . Whereas in the first two 
models 100 self interacting bodies (distributed between 0.5 and 4 
AU with small eccentricities) were integrated which then end up as 
planets, the last one starts with 400 embryos in an annulus between 
0.7 and 1 AU and - according to the authors_- repres ent mo re or 
less the outcome of the Grand Tack model dWalsh et al ll2012h . It is 
therefore appropriate to make such a study - which is orientated 
versus the collision of a Mars-sized object with the Earth - on the 
basis of these results. 
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Figure 1. Semi-major axes versus the mass of formed planets in an example Figure 2. Combined results of 16 different planet forming simulations; axes 

where the i nitial c onditions where chosen after the Grand Tack scenario like in Fig. |T| 

[Walsh et all2012l) . 


1.1 A possible formation of a ‘Theia’ in the early planetary 
system 

It is well known that the outcome of computations of the devel¬ 
opment of the early Solar System depends highly on the initial 
conditions after the gas in the disk disappeared. Several different 
approaches lead to different ‘planetary systems’ although all these 
attempts have b een undertaken to understand the architecture of 
our system (e.g. l lHansen 1 l2009l : llzidoro et al 112014) ). Our assump- 
tions is based on the outcome of the Grand Tack IWalsh et al l2012h 
where after the inward migration of Jupiter and Saturn a later out¬ 
ward migration triggered the formation of the terrestrial planets. 
In our numerical integrations we started with these two gas giants 
in their actual position and 100 planetesimals randomly distributed 
between 0.4 AU < flpianetesimai < 2.7 AU with masses in the order 
of the Moon. In Fig.jljwe plotted the results of one out of 16 sim¬ 
ulations where the architecture of this simulated planetary system 
turned out to be close to the one of our Solar System. Nevertheless 
there are two important differences: an additional planet (about the 
size of Mars) between the ‘Earth’ (here with only 60 % of its actual 
mass) and a ‘Mars' (with more than the double of its actual mass); 
the planet at 0.5 AU can be seen as a Venus equivalent. While this 
is just one example of resulting configurations we use it to moti¬ 
vate our choice of initial conditions described in the section below: 
a Mars-sized planet between Earth and Mars. The example chosen 
is not fully artificial when we look at Fig. [2] where all the simula¬ 
tions are combined into one graph which shows the accumulation 
of planets around 1 AU. 


2 DYNAMICAL MODEL AND NUMERICAL SETUP 

Our dynamical models for the early Solar System where chosen in 
the following way: 

• Model 1 [Ml], consisting only of the terrestrial planets ex¬ 
cluding Mercury and an additional planet in between the orbits of 
Earth and Venus. 

• Model 2 [M2], consisting of the terrestrial planets excluding 
Mercury, the gas giants Jupiter and Saturn, and an additional planet 
in between the orbits of Earth and Mars. 

Because our tests have shown that for the additional planet between 
Venus and Earth the influence of Jupiter and Saturn can be ne¬ 
glected we just used a 5-body problem for Ml. For both models 


we adopted as initial conditions the current osculating elements of 
the planets and added a Mars-sized object, Theia, with semi-major 
axes 0.8 AU < «xhcia < 0.94 AU (Ml) and 1.06 AU < an, c i a < 
1.4 AU (M2), respectively. Between scenarios we varied Theia’s 
semi-major axes with dczxheia = 0.005 AU. The reason of taking 
the orbital elements of the planets as they are today is the fol¬ 
lowing: we know that for the billio ns of years into the past the 
orbits were the same (e.g. lLaskarl Il996l) ) or only slightly differ¬ 
ent. The additional planet - the possible impactor - should then 
be in a quasi stable orbit after the formation of the planetary sys¬ 
tem. To find such an orbit we did not vary the eccentricity nor 
the inclination and set them to e m \ = 0.075 and i = 2°. For the 
other orbital elements of Theia we used randomly chosen values 
between 0° and 360°. For every fixed semi-major axis 25 such os¬ 
culating elements were computed as initial conditions. For achiev¬ 
ing the highest possible precision with respect to the collision an¬ 
gle and velocity the integration method was the one we hav e used 
for many years for most of our computations ( Dvora lTet al II2003L 
120151 : iGaliazzo et al ll20L3h , The Lie-integration has an automatic 
step-size control and is well adapted for such kind of computations 
jHanslmeiem &Dvorak II 1 984] : lEggl & Dvorakll2oTcl : iLichtenegger I 


ll984l:lDelva|ll985t) . With regard to the formation process of the 

terrestrial planets and the estimated time of co llision of Theia with 
the Earth of 95 ± 32 Myr ijacobson & Morbidelli I li2014l) ) the inte¬ 
gration time was up to 50 Myrs. 


3 EVOLUTION OF TWO SELECTED SAMPLES 

To demonstrate the sensitivity of the dynamical evolution in this re¬ 
gions we show a stable seeming orbit and another one very rapidly 
becoming unstable. Both had initial conditions outside the Earth or¬ 
bit with only slightly different semi-major axis, and the same initial 
conditions for the other orbital elements. 

In Fig. [3] and Fig. [5] we show two examples of the develop¬ 
ment of the semi-major axes of a fictitious Theia and the terrestrial 
planets, one for a stable orbit and another one which is chaotic very 
soon. The very regular variation of the eccentricity of this stable or¬ 
bit of Theia is depicted in Fig.[4](0 < e-rheia < 0.08) which ends in a 
sudden increase up to e = 0.15, a close approach and even a colli¬ 
sion with the Earth after 17 Myrs. In the other example (Fig. 0 we 
show a chaotic orbit suffering from multiple close encounters with 
the Earth visible through the chaotic signal of Theia’s semi-major 
axis already after 1 Myr. Caused by a sequence of very close en¬ 
counters between Earth and Theia the orbit of the fictitious planet 
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time in 1000 years 

Figure 3. Time development of the semi-major axes of the terrestrial planets 
and additional planet Theia leading to a collision with the Earth. 



Figure 6. Time development of the eccentricity (y-axis) of a highly chaotic 
orbit of Theia from Fig.[5]with a subsequent escape after a close encounter 
with the Earth; the eccentricity of Mars is the regular curve close to e = 0.1 



Figure 4. Time development of the eccentricities of the Earth, Venus, and 
in addition the one of Theia (in red) from the example in Fig. [5] Please see 
text for more. 



Figure 7. Schematic view of the collision parameters impact angle and im¬ 
pact velocity for an Earth-colliding small planet. 


jumps inside the orbit of the Earth after 2.2 Myrs. Also the Earth is 
suffering from these repeatedly close encounters and moves a lit¬ 
tle outside which can be explained by the conservation of the mo¬ 
mentum of this planet pair. After a capture close to the 6:5 MMR 
between them (2.2 Mrys) a final collision with the Earth ends the 
lifetime of Theia. Mars and Venus are not affected at all in this ex¬ 
ample. 



Figure 5. Time development of the semi-major axes of a highly chaotic 
orbit leading to an escape after 3 Myrs. 


4 IMPACT PARAMETERS 

The goal of these investigations is to determine the collision prob¬ 
ability for an additional Mars-sized planet Theia which could have 
formed during the early formation (e.g. Grand Tack scenario) be¬ 
tween the Earth and Mars or between Earth and Venus. In addition 
to this we study the impact velocity as well as the impact angle 
(Fig- 0 - expressed by the impact parameter in diff erent studies 
dLeinhardt & Stewartll2012l : IMaindl & Dvorak ll 201 4|) . 

Together with the mass and composition of the two colliding 
bodies these two quantities determine the outcome of such a cosmic 
catastrophic event. In Fig. [8] we show the different impact parame¬ 
ters according to our simulations, where one can see that most of 
the encounter velocities are just above the escape velocity of the 
two planets (v esc ~ 9km/s) and only one impact happens with a 
15 percent higher velocity. The distribution of the impact angles is 
asymmetric with slightly more head-on collisions; we can explain 
it by the fact that close to the Earth Theia suffers from a strong ac¬ 
celeration toward the center of the larger planet (we remember the 
mass ratio /j = 0.1). 

The actual collision outcome in terms of fragmentation 
strongly de pends on the impact velocities, impact angles, a nd the 
mass ratio dLeinhardt & Stewart 11 201 2l ; IMaindl et alll2014h . The 
number of fragments varies with the impact angle such that even for 
relatively high impact velocities strongly inclined collisions may 
lead to two major survivors (hit-and-run) whereas head-on impacts 
may destroy the involved bodies. In our scenarios the low impact 
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angle [deg] 

Figure 8. Collision parameters impact angle versus impact velocity for the 
Earth-colliding Theias. 



aiheialAU] 

Figure 10. Schematic view of the domain of the terrestrial planets with 
respect to the semi-major axis. In the regions b and d we plotted the number 
of observed collisions with the Earth; for more see text. 
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Figure 9. Initial condition diagram: differences in mean logitudes with re¬ 
spect to the Earth versus the semi-major axis. The color points are coded 
with respect to their time of stability. 


«Theia 

CT-mean 

CT-max 

CT-min 

n 

0.875 

3.41 

6.648 

1.555 

3 

0.880 

6.13 

6.129 

6.129 

1 

0.885 

12.87 

18.542 

2.892 

3 

0.890 

8.52 

12.691 

4.342 

2 

0.895 

7.02 

13.973 

3.332 

4 

0.900 

8.21 

23.669 

0.313 

3 

0.905 

1.83 

4.528 

0.114 

3 

0.910 

2.06 

2.904 

1.216 

2 

0.915 

1.92 

3.747 

0.090 

2 

0.920 

4.22 

6.168 

1.628 

4 

0.925 

9.76 

15.685 

0.057 

3 

0.930 

0.26 

0.626 

0.051 

3 

0.935 

0.02 

0.024 

0.024 

1 

0.940 

0.97 

4.872 

0.002 

9 


Table 1. Collisions in region b, see text for details. 


velocities suggest surviving bodies and probably a Moon forming 
from a debris ring. 


5 THE OVERALL PICTURE 

For the integrated several thousands of orbits and the chosen grid of 
semi-major axes a fixed eccentricity and a fixed value of the incli¬ 
nation was chosen; the initial mean longitude was varied randomly 
(see Sect. [2j. In Fig.[9]we show the combined results for all com¬ 
puted orbits by plotting the initial difference in mean longitudes 
(rfrheia - d Earth ) versus semi-major axes and color coding the ‘sta¬ 
bility character’ of the orbit which was classified as follows: red to 
light blue circles according to the escape time. These color points 
stand for an unstable orbit either due a close encounter with Venus 
(inner part) or due to a close encounter with Mars or even with the 
Earth. We note such an orbit as one typical example is depicted in 
Kg. m Dark blue circles stand for stable orbits for the whole inte¬ 
gration time; the black dots stand for a collisions with the Earth. 

According to the chosen initial conditions in semi-major axis 
we divided the whole domain where the terrestrial planets move in 
5 regions: (a) aTheia < 0.875 AU, (b) 0.875 AU < aTheia < 0.94 AU, 
(c) the region around the Earth, (d) 1.06 AU < a Theia < 1.165 AU, 
and (e) aTheia > 1.165 AU. In regions b and d we plotted the num¬ 
ber of collisions with the Earth for each of the 25 initial mean lon¬ 


gitudes per initial aTheia (Fig. m. We observe that the maximum 
number of collisions with the Earth is about 36 % (9 out of 25 ini¬ 
tial conditions), which occurs once for regions b (a Theia = 0.95 AU) 
and d (fl-rheia = 1.075 AU), respectively. 

In regions a and e many mean motion resonances (MMR) with 
the Earth act to destabilize the orbits of a hypothetical Theia, see 
Fig.[TT]which also shows the most important MMRs together with 
the longest duration of stability of the collision orbits in regions b 
and d. One can see that the number of Earth collisions increases 
with the decreasing distance to the Earth. We did not make compu¬ 
tations close to the Earth (inside and outside its orbit) - region c - 
with the exceptions of some sample orbits, which have shown that 
almost all are destabilized very soon after a close encounter or even 
a collision. Because our interest was to find collisional orbits which 
were stable for several tens of millions of years in between Earth 
and Venus and also between Earth and Mars we do not show these 
results. It is easy to understand that only then an additional planet 
may have formed in the early formation stages of our planetary 
system; we showed even an example from our own computations 
of the formation of the planets in FigQ] where a Theia like planet 
was formed outside the Earths orbit. 

It is obvious that close to the inner region c the stability time is 
short because the Earth is relatively close by, whereas we find times 
up to 30 Myrs of stability for an orbit before an escape in region d. 
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Figure 11. Mean Motion Resonances of Theia with the Earth in regions b 
and d. Also the longest time before a collision versus semi-major axis is 
plotted here; for more see text. 



a Th.ia I AU 1 


Figure 12. Mean time to collisions in Myrs versus semi-major axis of Theia 
in AU. 


The results for the mean time before a collision occurs are shown 
in Fig. |T2] It is obvious, that close to the inner region c this time 
is short because the Earth is relatively close by, whereas this mean 
collision time increases in regions b and d. 


6 CONCLUSIONS 

Thousands of orbits of a hypothetical additional planet Theia in 
the early phase of the Solar System were integrated and classified 
according their dynamics. The aim was to find an orbit stable for 
sufficiently long time (tens of Myrs) inside or outside the orbit of 
the Earth which then could lead to a collision building a companion 
of the Earth. The impactor was assumed having the mass of Mars 
which then - after a collision - could lead to an additional body 
like our Moon. 

Tabs.|7]and[2]summarize our results with respect to the mean, 
maximum and minimum collision times for the inner region b and 
the outer region d, respectively. The last columns show the num¬ 
ber of Earth collisions of Theia, which is always smaller than 40 
percent of the 25 integrated orbits for one fixed semi-major axis. 
Finally, Tab. [3] shows the overall statistics. Approximately 50 % 
of all integrated orbits (in total about 2000, the different tests in 
regions a and e included) turned out to be stable up to the cho¬ 
sen integration time of up to 50 Myrs. The same almost 50% in 


^Theia 

CT-mean 

CT-max 

CT-min 

n 

1.065 

0.17 

0.235 

0.013 

2 

1.070 

0.52 

1.255 

0.008 

4 

1.075 

0.86 

4.518 

0.010 

9 

1.080 

3.60 

10.471 

0.017 

8 

1.085 

11.10 

0.314 

0.104 

4 

1.090 

3.22 

18.218 

0.023 

7 

1.095 

8.21 

14.190 

0.205 

6 

1.100 

15.20 

8.838 

0.117 

4 

1.105 

2.30 

17.066 

0.426 

4 

1.110 

2.62 

6.225 

0.392 

3 

1.115 

2.59 

3.580 

1.259 

3 

1.120 

5.28 

15.63 

1.529 

5 

1.125 

23.43 

23.435 

23.435 

1 

1.130 

- 

- 

- 

- 

1.135 

10.19 

14.778 

1.299 

3 

1.140 

10.12 

9.691 

5.27 

2 

1.145 

7.97 

18.064 

3.525 

7 

1.150 

- 

- 

- 

- 

1.155 

9.44 

9.443 

9.443 

1 

1.160 

31.19 

31.193 

31.193 

1 

1.165 

20.26 

30.895 

9.617 

2 


Table 2. Collisions in region d 


region 

flTheia in AU 

stable 

eject 

collision 

a 

0.750-0.875 

50.50 

47.25 

2.25 

b 

0.875-0.940 

47.08 

40.61 

12.31 

c 

0.940-1.060 

- 

- 

- 

d 

1.060-1.165 

44.73 

29.27 

26.00 

e 

1.165-1.350 

52.16 

45.12 

2.72 


Table 3. Statistics of collisions (units percents) 


these two regions escape due to close encounters either with Venus 
or Mars (we did not count the number of Venus and Mars colli¬ 
sions) and only 2% suffered from impacts on Earth. In regions 
b again approx. 50 % turned out to be stable, but the number of 
Earth colliders was increased to 12 %. Region d is the best candi¬ 
date for a collision of Theia with the Earth after millions of years 
of stability: 26 % of such colliders were observed here and espe¬ 
cially around 1.16 AU (compare Fig.l 1 1 1 the chance is high that the 
Moon-producing planet Theia was formed here together with the 
other terrestrial planets (compare Fig.Q}. We plan to make more 
of these numerical experiments on a finer grid and for longer inte¬ 
gration times and furthermore we will combine these results with 
detailed computations of the collisions (SPH codes); all this will 
bring us a step further in the knowledge of the formation of our 
Moon. 


7 ACKNOWLEDGEMENTS 

Special thanks for financial support are due to the FWF, project 
S-11603-N16. 


REFERENCES 

Asphaug, E Annual Rev.of Earth and Planetary Sciences, 42, 551 
(2014) 


© 2015 RAS, MNRAS 000,[Q-?? 










































6 R. Dvorak 


Canup, R., Asphaug, E., Nature, 41, 6848, 708 (2001) 

Canup, R„ Icarus, 168, 433 (2004) 

Canup, R. M„ Icarus , 196, 218(2008) 

Canup, R„ Nature, 504, 27 (2013) 

Cameron, G.W ., Icarus, 126, 126 (1997) 

Cameron and Ward, Lunar Planet. Sci. Conf. (1976) 

Cuk, Matija; Stewart, Sarah T. Science, 338, 1047 (2012) 

Delva, M„ A&rA 60, 211 (1985) 

Dvorak, R., Pilat-Lohinger, E., Funk, B., Freistetter. F., ASrA 410, 
L13 (2003) 

Dvorak, R., Maindl, T. I., Siili, A., Schafer, C., Speith, R., Burger, 
C„ 2015pees.confP...ID (2015) 

Eggl, S. and Dvorak, R., Lecture Notes in Physics, 790, 
431(2010). 

Galiazzo, M. A„ Bazso, A„ Dvorak, R..POSSX4, 5 (2013) 
Hansen, B. M. S„ ApJ 703, 1131 (2009) 

Hanslmeier, A.; Dvorak, R. AbA 132, 203,(1984) 

Hartmann, W.K., Davis, D.R., Icarus, 24, 504 (1975) 

Izidoro, A., Haghighipour, N., Winter, O. C., Tsuchida, M., ApJ, 
782, 31 (2014) 

Jacobson, S.A.,Morbidelli, A.,Raymond, S.N., O’Brian, D.P., 
Walsh, K.J. Rubie, D.C., Nature, 508, 84 (2014) 

Jacobson, S.A.,Morbidelli, Phil.Trans.R.Soc.A. 372, 0174 (2014) 
Kaib, N.A., Cowan, N.B., Icarus, 252, 161 (2015) 

Laskar, J„ CMDA 64, 115 (1996) 

Eeinhardt, Z.M. & Stewart, S.T., ApJ, 745, 79 (2012) 
Eichtenegger, H„ CMDA 34, 357 (1984) 

Maindl, T.I. & Dvorak, R. (2014), IAU Symposium 299, 370 
(2014) 

Maindl, T.I., Dvorak, R., Schafer, C., Speith, R., IAU Symposium 
310, 138 (2014) 

Nakajima, M„ & Stevenson, D. J„ arXiv:1506.04853 (2015) 
Quarles, B.L., Lissauer, J.J., Icarus, 248, 318 (2015) 

Reufer, A.; Meier, M.M. M., Benz, W., Wieler, R., Icarus, 221, 
296 (2012 

Walsh, K. J., Morbidelli, A.. Raymond, S. N. O’Brien, D. P, Man- 
dell, A. M„ 2011, Nature, 475, 206 (2012) 


© 2015 RAS, MNRAS 000,[Q-?? 


